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ABSTRACT 

Large nebulae (>50 kpc) emitting strongly in Lya (also known as Lya "blobs") are likely signposts 
of ongoing massive galaxy formation. The relative rarity of these sources and their discovery in well- 
studied galaxy overdensities suggest that they may be associated with regions of high galaxy density. 
One of the largest Lya nebulae, discovered at a redshift of z w 2.7 via its strong mid-infrared emission, 
provides an unbiased test of this association. We have carried out a deep intermediate-band imaging 
survey for Lya-emitting galaxies (LAEs) within a 30'x26' field of view around this Lya nebula. 
This is the first study of the environment of a Lya nebula found without a priori knowledge of its 
surroundings. We find that the nebula is located in an overdense region, at least 20x50 hf comoving 
Mpc in size, showing a factor of ^3 LAE number density enhancement relative to the edge of the field. 
Given the predicted number of such overdensities, we rule out the possibility of a chance coincidence 
at the <1% level. This study, in conjunction with previous work, provides strong confirmation of the 
association between the largest Lya nebulae and overdense regions of the Universe. 
Subject headings: galaxies: formation — galaxies: high-redshift — large-scale structure of universe 



1. INTRODUCTION 

Studies of massive galaxy populations show that the 
most massive galaxies are in place and in posses- 
sion o f the majority of the i r stellar mass by z ~ 

1 - 2 ([McCarthy et alj 12004 Ivan Dokkum et al l 12004 



iDaddi et alJl200alBundv et ail2005tlBrown et aljl2007f ). 
Thus, while dark matter halos in a ACDM cosmol- 
ogy build up hierarchically with the gradual accretion 
and merging of smaller halos, the most massive galax- 
ies within that context likely have more dramatic ori- 
gins. The details of this process are uncertain, and 
ideally we would like to study sites of ongoing massive 
galaxy formation. Lya nebulae or Lya "blobs" (LABs) 
- large (>50 kpc) clouds of gas emitting strongly in 
Lya (~ 10 44 erg s _1 ) - provide such an opportunity. 
The large Lya equivalent widths and the association 
with galaxy populations such as Lyman break galax- 
ies (LBGs) and submillimcter galaxies (SMGs) strongly 
suggest that LABs are sites of ongoing galaxy forma- 
tion. They have been found in small numbers (only 16 
>50 kpc LABs arc known) around z ~ 2-3, a key epoch 
of black hole and galaxy growth (ISteidel et al 
hereafter S00; iFrancis et al.l 120011: iPalunas et al 



200 



200 



Matsuda et "ail 2004: D ev et al.ll20 05: Nil sson et alJl2006t 
Smith h Jarvisll2007t iGreve et al.ll2007fh LABs span a 
range in size and surface brightness from the scales of 
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Lya-emitting galaxies (LAEs) to the large st LABs known 
(~15 kpc, ~27 mag arcsec -2 ; e.g., iMatsuda et al.l 
I2004T) . LABs are similar to the large Lya halos ob- 
served in the ove rdensities around higher redshift ra- 
dio g alaxies fe.g.. IVenemans et al.l 120071 : lOverzier et al.l 
2008) but are found in radio quiet environments. The 
dominant power source in LABs is difficult to deter- 
mine; among the limited sample of large LABs, there 
is evidence for embedded AGN, starburst-driven super- 
winds, gravitational cooling radiation, and spatially- 
extended star formation, all of which may play a role 
in powering the Ly a emission (Prescott et al. 2008, 
in pr e paration; S00: IChapman et al.ll2004 iNilsson et al.l 
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Taniguchi h Shioval l20qfljMatsuda et all 12004 
2005; Matsuda ct al. 2007). 



Most known LABs, including two of the largest cases, 
have been discovered via narrow-band sur veys, often by 
targeting known galaxy overdensities fS00:|Francis et al.l 
1200 It IPalunas et al.ll2004 IMatsuda et al.ll2004 ). Follow- 
up narrow-band imaging of the SSA22 region surround- 
ing the S00 LABs revealed fainter LABs associated with 
the sa me overdensity traced by the LAEs (Ma tsuda et al.l 
|2004[ ). This suggests that LABs may be confined to over- 
dense regions, as would be expected for site s of massive 
galaxy formation. A blank field survey by ISaito et"aTI 
(2006) supported this claim, finding no LABs with sizes 
greater than ~30 kpc and an order of magnitude lower 
number density of LABs relative to that found within 
the SSA22 galaxy overdensity. However, the association 
between LABs and galaxy overdensities may be mislead- 
ing, as a truely systematic wide-area search has yet to 
be completed. A thorough environmental study has only 
been done for the S00 LABs, which were found by tar- 
geting a known galaxy overdensity. 

In contrast, one of the largest LABs uncovered recently 
was found by entirely different means. While conducting 
a study of mid-infrared ( 24/nm) sou r ces de tected by the 
Spitzer Space Telescope, iDev et all (|2005h discovered a 
LAB at z = 2.656 within the NOAO Deep Wide-Field 
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Survey Bootes field (NDWFS 7 ; Uannuzi fc DevlH999l) . 
Follow-up observations revealed complexity typical of 
the LAB class: the region hosts a buried AGN, many 
young galaxies, an LBG, and diffuse Hell and contin- 
uum emission suggestive of spatially-extended star for- 
mation (Prescott et al. 2008, in preparation). One of 
the largest and most luminous LABs known, and found 
without any a priori knowledge of its surroundings, this 
source represents a unique opportunity to perform an un- 
biased, complementary test of the association between 
LABs and overdense regions. 

In this work, we present the first results from an ongo- 
ing deep intermedia te-band Ly c t imag ing survey of the 
environment of the iDev et all (|2005 ) LAB (hereafter, 
LABd05) and report on the spatial distribution of the 
LAEs in the immediate vicinity. A detailed analysis of 
the multiwavelength properties of LABd05 and the prop- 
erties and clustering of the LAEs in the region will be pre- 
sented by Prescott et al. (2008, in preparation). We as- 
sume the standard ACDM cosmology (S1m=0.3, Oa=0.7, 
h=0.7); the angular scale at z = 2.656 is 7.96 kpc/". All 
magnitudes are in the AB system. 

2. OBSERVATIONS & REDUCTIONS 

We obtained deep imaging of the field around LABd05 
using the Subaru tele scope and the SuprimeCam imager 
(jMivazaki et al.ll2002f ) on U.T. 2007 May 10-14 and June 
17. The survey covers 0.22 deg 2 in an intermediate-band 
filter, IA445 (A c «4458A, AX FWHM »20lA), centered 
on the Lya line at the redshift of the nebula; this cor- 
responds to a comoving volume of 4.27xl0 5 h^ 3 Mpc 3 
(52x45x180 hfg 1 Mpc). Conditions during the May ob- 
servations were variable (clouds, variable seeing 0.7"- 
1.2") and good in June (clear with 0.7" seeing). We 
obtained a total of 3 hours of observations. 

We reduced the data using the SDFRED software 
ijYagi et all |2003 lOuchi et all l200l . The dat a were 
overscan-subtracted and corrected for geometric and at- 
mospheric distortions. We generated the sky flat using 
object frames that were free of bright stars in combina- 
tion with other images taken in the same intermediate- 
band filter (Y. Taniguchi, private communication). Small 
portions of the SuprimeCam field of view are vignetted 
by the autoguider probe; the affected areas were masked, 
as were bad columns, bright star ghosts, and satellite 
trails. Images were aligned and scaled using common 
stars and then combined using a clipped mean algorithm, 
which successfully removed cosmic rays. Of the Suprime- 
Cam field-of-view (0.26 deg 2 ), 73% was usable (the re- 
mainder being of lower signal-to-noise along the edge of 
the field or in the vicinity of bright stars). 

The limiting magnitude of the stacked image is 28.3 AB 
mag (lcr, 2" diameter aperture), calculated using 10,000 
random apertures. An approximate magnitude zeropoint 
was calculated from observations of the standard stars 
BD+25d4655 and Feige 34. For the NDWFS broad-band 
imaging, the limiting magnitudes are Bwiimit—^7.9 mag, 
Rumtt=27-l mag, / KmW =26.0 mag (lcr, 2" diameter aper- 
ture). 

A portion of the field was observed in a custom U- 
band filter (U d ; A «3590A, A\ FWHM -116A) using the 

7 This research draws upon data provided by Dr. Buell Jannuzi 
and Dr. Arjun Dey as distributed by the NOAO Science Archive. 



Mayall 4m Telescope over 6 nights (U.T. 2007 June 8- 
13). These data will be described elsewhere (Prescott et 
al. 2008, in preparation) but here provide a useful check 
on interlopers in our z w 2.7 LAE sample. 

3. CANDIDATE SELECTION 

We used Source Extractor (jBertin fc Arnoutill996T ) to 
select a sample of «38,600 sources detected in the 7A445 
band down to the ha limit of 7.4445=26.5 mag (2" di- 
ameter aperture; Lia445{z — 2.656) = 1.5xl0 42 erg s _1 ) 
with the following search parameters: at least 5 contigu- 
ous pixels, a threshold of 2a per pixel, and a Gaussian 
filter matched to the seeing (FWHMw0.8 arcsec). We 
measured matched aperture photometry using 2" diam- 
eter apertures from the 771445, Bw, R, and 7 imag- 
ing datasets, which were registered and resampled to 
match the Bw astrometry and pixel scale (0.258 arc- 
sec/pix). Aperture corrections were neglected (ral.08 
for an unresolved source). Line-emitting sources are 
strongly detected in IA445 relative to Bw, i.e., they 
have large negative IAAAh—Bw colors relative to the 
normal galaxy locus (see Figure [l}. We removed bright 
stars (7A445< 25.0) using the CLASS_STAR param- 
eter in Source Extractor (>0.91) and employed a cut 
of IAA4:5~B W < -0.85 mag yielding 1500 candidates. 
Shifting the I AAAh—Bw cut by ±0.2 mag causes no sig- 
nificant change to the main results presented in Section 4. 
For a F v oc v continuum source, this corresponds to an 
observed equivalent width (EW) cut of W bs >148 A. 
LAEs are kno wn to be young, with estimated ages 
of 4-200 Myr (iFinkelstein et al. 1 120071 lLai et alJ 120071 : 
iGawiser et al.l I2007D . For the case of a young galaxy 
(25 Myr old si mple stellar population, solar metallicity, 
Chab rier IMF; iBruzual fc Chariot! 120031 : iTremonti et"atl 
l2004f) at z = 2.7 with standard intergalactic absorption 
( Madaulll995l ), this is equivalent to a rest-frame EW cut 

Of W rest >50A. 

We expect our LAE candidate sample to be con- 
taminated by high EW, low-redshift [OII]AA3727,3729- 
emitting galaxies at z ks 0.2. Since the 771445 bandpass 
lies on the red side of the Bw filter (Figure Q] see in- 
set), we also expect contamination from higher redshift 
galaxies (z > 2.9 — 4.0) for which the Lyman limit has 
entered the Bw filter, thus depressing the Bw flux rel- 
ative to the 7A445 flux. Using the pu blicly available 
NDWFS imaging (| Jannuzi fc Devi [l999t ). we employ a 
cut of Bw~R <0.8 to remove both contaminant pop- 
ulations (see Figure [2]). The final sample of 785 LAE 
candidates corresponds to a mean LAE surface density 
of «4200 deg~ 2 . The properties and sizes of the LAE 
sample will be discussed in an upcoming paper (Prescott 
et al. 2008, in preparation). There are no other large 
(>50 kpc) LABs in the vicinity of LABd05. 

Sources at these redshifts should show very little 
flux in the t/d-band, which straddles the Lyman limit 
{X rest «970-1000A at z « 2.656). Using the 25 Myr 
old model above to represent the typical LAE con- 
tinuu m shape, th e z ~ 3 galaxy luminosity function 
from lReddv et al.l (l2007h . and standard intergalactic ab- 
sorption (|Madaul Il995f ). we predict that ^1-2% of the 
LAE sample should have Ud detected at the 5cr level 
(f/ d =25.3). Only 1 LAE candidate (<1%) is detected 
with Ud <25.3, giving us confidence that we have selected 
a clean sample of LAEs. More sophisticated interloper 
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rejection would require deeper Ud imaging over the entire 
field. 

4. RESULTS AND DISCUSSION 

Figure [3] shows the spatial distribution of LAEs in the 
vicinity of LABd05. We find an overabundance of LAEs 
in the immediate vicinity of the LAB and clear evidence 
of an elongated overdense structure, with what appears 
to be a ~10 h^ 1 Mpc underdensity at [218.36,33.07]. In 
Figure HI we plot the surface density of LAEs versus dis- 
tance from LABd05. This shows a peak overdensity of a 
factor of 3.0 relative to the edge of the field, or a factor 
of 2.6 when averaged over the central 10 h^ 1 Mpc ra- 
dius. For comparison, the dashed line shows the expected 
LAE surface de nsity estimated f r om th e z ~ 3.1 luminos- 
ity fu nctions of iGronwall et al.1 (|2007l ) and lOuchi et al.1 
(l2007h assuming a uniform rcdshift distribution across 
our bandpass and no evolution. The overdensity spans 
at least «17x47 hyp 1 Mpc (comoving). A study of the 
LAE population around the two S00 LABs, which are of 
comparable size and luminosity to LABd05, also found 
a factor of ~3 o verdensity relative to t he field that was 
~60 Mpc across (jHavashino et al.ll2004h . Thus, although 
it was found without any a priori knowledge of its sur- 
roundings, LABd05 appears to reside in a similarly over- 
dense environment. 

Without follow-up spectroscopy, we cannot know the 
true redshift distribution of the sample within the 
intermediate-band filter. Assuming a uniform distribu- 
tion across the filter (Az « 0.17, 180 hyp 1 Mpc, comov- 
ing, along the line-of-sight), our survey yields a mean 
number density of p «2.1xl0~ 3 h^ 3 Mpc~ 3 for LAEs 
with L/A445 >1.5xl0 42 erg s -1 cm~ 2 . The inner con- 
tour in Figure [3] corresponds to p «2.8xl0 -3 hy 3 Mpc -3 
while at the edge of the field p wl.2xl0 -3 h^ 3 Mpc -3 . 
Given the presence of an overdensity, we expect the true 
redshift distribution to be significantly narrower. If we 
assume the overdensity resembles that hosting the S00 
LABs, which extends 40 hyp 1 Mpc co moving (22% of 
the fi lter width) along the line-of-sight (jMatsuda et al.l 
2005), and that it is centered within the IAAA5 filter, 
the inner contour corresponds to a number density of 
p «12xl0 -3 hy 3 Mpc" 3 . Alternately, if the structure is 
cylindrical in shape with a radius of 10 hy Mpc, the 
corresponding number density at the inner contour is 
p «25xl0~ 3 hfo 3 Mpc" 3 . 

We estimate the expected frequency o f such over- 
dense structures using the galaxy catalog of iBower et al.1 
(2006) , which is ba sed on the Millennium Simulation 
(jSpringel et al.|[2005l ). We approximate our observational 
set-up by sampling the simulation volume randomly with 
cylinders that are 10 hyp 1 Mpc in radius (the size of 
the overdensity peak) and 180 hyQ 1 Mpc in depth (the 
full sp an of the filter). We assumed a Lya /Ha ratio 
of 10 (jOsterbrockl 11989] ) and scaled from the predicted 
Ha luminosity to select model LAE galaxies down to the 
Lya limit of our observations (under the assumption that 
Lya dominates the measured IAAAh flux) at a redshift of 
z w 2.7. The number density of model LAEs in this case 
is 10.8 xlO -3 hy 3 Mpc -3 , which is a factor of ~5 higher 
than the mean density we observe. Taken at face value, 
this could imply that LAEs have a duty cycle of —20% or 



that they are a younger, less massive subset of this pop- 
ulation. If we restrict the model galaxies to be younger 
than the sample's median age (<172 Myr) and less mas- 
sive than the median stellar mass (<9xl0 8 M Q ), the 
number density of model LAEs is 3.5 x 10~ 3 hy- 3 Mpc" 3 . 
In either case, we uncover overdensities of greater than 
a factor of 2 at a rate of ssO.3%. Therefore, within the 
large span of our filter, which will tend to average out 
inhomogeneities and reduce the signal, our imaging sur- 
vey had only a ~0.3% chance of randomly uncovering 
such an overdense region if LABs and overdensities are 
independent phenomena. The space density of large Lya 
nebulae is very uncertain, but at the high end is the range 
quoted for the S0 Lya nebulae; -3-400 xl0~ 6 Mpc -3 
(jSaito et al.ll2006l) . Taking these values (equivalent to a 
—17-100% chance of finding one large LAB within the 
peak of the overdensity) , the likelihood of a chance coin- 
cidence between a factor of >2 overdensity and a large 
Lya nebulae would be <0.05-0.3%. Preliminary results 
from more recent systematic Lya nebulae surveys hint 
that their true space density may be orders of magni- 
tude lower (e.g., Prescott et al. 2008, in preparation), 
making the likelihood of chance coincidence vanishingly 
small. 

In terms of the relevant emission mechanisms the Lya 
nebula class appears to be a highly heterogeneous mix, 
and this diversity could in principle derive from environ- 
mental differences. The largest Lya nebulae (>100 kpc), 
including the case studied here, often show evidence 
for obscured AGN and extended star format ion (e.g., 
iMatsuda et"alll2007t iBasu-Zvch fc Scharj[200l Prescott 
et al. 2008, in preparation) and have received the most 
scrutiny in terms of their properties and environments. 
Presumably, the somewhat smaller "coolin g" LABs (e.g., 
iNilsson et all 120061 : ISmith k. Jarvisl l2007f ) must also re- 
side in dense regions with sufficient gas supply, but thor- 
ough environmental studies of these sources have, to our 
knowledge, not yet been completed. 

5. CONCLUSIONS 

The discovery of a large Lya nebula at z « 2.7 via 
its strong mid-infrared emission has provided an unbi- 
ased test of the association between these rare sources 
and galaxy overdensities. Using deep Lya imaging of 
the environment surrounding this LAB, we identify 785 
LAE candidates and find evidence for a factor of —3 LAE 
overdensity which spans 20x50 Mpc (comoving). This is 
comparable to what is found in the vicinity of the well- 
known S00 Lya nebulae. We rule out a chance coinci- 
dence at the <1% level. In conjunction with previous 
work, these results point conclusively to a strong asso- 
ciation between the largest Lya nebulae and overdense 
regions of the Universe. 

We are grateful to Masafumi Yagi and Yoshiaki 
Taniguchi for data reduction advice and assistance. We 
thank Naveen Reddy for the use of data from an ongo- 
ing Ud imaging program and the referee for constructive 
comments on the manuscript. M. P. was supported by 
an NSF Graduate Research Fellowship. Support for this 
work was provided by NASA (HST-GO10591) and the 
National Science Foundation (Grant No. 0714311). 
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Fig. 1. — IA4A5—Bw color-magnitude diagram for sources detected in the 7A445 image (dots). Initial candidates were selected to meet 
the color cut I A445— Byy < —0.85 (dashed line, which corresponds to an observed equivalent width of W t, 3 >148A) and to be brighter 
than the 5<r limiting magnitude in IA445 band (//14455 CTi Hmit=26.6; dotted line). LABd05 is shown as a star, and non-detections in the 
Bw band are set to the Bw limiting magnitude. Typical error bars are shown. The inset shows the effective transmission curves for the 
Bw and 7/1445 filters and the wavelength of Lya at the redshift of the system (dotted line). 
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Fig. 2. — Bw—R vs. R — I color-color diagram of the 1500 line-emitting candidates (ope n and filled circles). A typical error bar 
is shown. Expected color tracks for low redshift galaxy templates are shown for O.K z <0.3 IILeitherer et al.lll996T) ; the color track for 
high redshift LAE/LBG c ontaminants is based on a young galaxy template at 2.9< z <4.0 (25 Myr old simple stellar population, solar 
metallicity, Chabrier IMF; Bruzual & Chariot 2003) with standard intcrgalactic absorption (Madau 1995). To remove both low and high 
redshift interlopers, candidates were required to have B-^—R <0.8. The final sample of 785 LAE candidates is shown (filled circles). 
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Fig. 3. — Spatial distribution of LAE candidates (filled circles). After smoothing on 3.2'x2.8' scales, isodensity contours were laid down 
at 0.7, 1.0, 1.8, and 2.3 times the average LAE density at the edge of the field. The four circle sizes represent four bins in IA445 luminosity: 
£/A445 < 2.5 x 10 42 , 2.5 x 10 42 < L IA445 < 5 x 10 42 , 5.0 x 10 42 < L IA445 < 1.0 x 10 43 , L IA445 > 1.0 x 10 43 . The position of LABd05 
is indicated with a star. Regions masked along the edge and due to bright stars are shown in grey. There is evidence for an extended 
structure stretching from SE to NW and a ~10 hfp Mpc void at [218.36,33.07]. 
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Fig. 4. — The radial density profile measured in circular annuli centered on LABd05. The LAB lies in an overdense region that is a factor 
of ~3 times the density at the edge of the fiel d of view. The dashed line represent th e predicted number density at z ss 3 if we assume a 
uniform redshift distribution and no evolution (Gronwall ct al. 2007; IOuchi et al.ll2007h . The shaded areas correspond to the range covered 
by lcr error b ars; all three param eters (a,L* ,d>*) from l^rorrw^IFetaLl 120071) were allowed to vary in turn (light shading), while in the 
comparison to lOuchi et al.l 120071) the faint-end slope a is fixed at —1.5 (dark shading). The predicted LAE number density is consistent 
with what we measured at the edge of the field. Since the redshift distribution of our sample can only be narrower than assumed, the 
overdensity measured within the structure is likely a lower limit. 



